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Abstract 



In this note, we study the eigenvectors and the scalar products the integrable 
long-range deformation of the XXX spin chain defined in [1]. The model is solved 
r~| ' exactly by algebraic Bethe ansatz, and it coincides in the bulk with the Inozemtsev 

spin chain. At the closing point it contains a defect which effectively removes the 
' wrapping interactions. Here we concentrate on determining the defect term for 

^ ■ the first non-trivial order in perturbation in the deformation parameter and how 

it affects the Bethe ansatz equations. Our study is motivated by the relation with 
^f) ' the dilatation operator of the Af = 4 gauge theory in the su{2) sector. 

rn ' 

I Introduction: Long range spin chains. In this paper we consider long- 

■ range integrable deformations of the XXX spin-1/2 spin chain. There are several 

methods to turn a nearest-neighbor spin chain into a long-range ones. One method, 
used in [2, 3], is to deform the conserved charges according to 



X: ^Qr{X)=i[xW,Qr{X)] , (1) 



where A the deformation parameter and xW some deformation operator, which in 
the simplest case can be the boost operator of one of the higher conserved charges. 
The integrated version of the formula above gives 

Qr{X) = S{X)Qr{0)S-HX) , with ^S-^=ix{X). (2) 

dX 

This procedure can be realized on chains of infinite length, but on chains of finite 
size there are two complications. First, the definition of the operator S{X) might 
not be compatible with periodic boundary conditions, so the resulting chain will not 
be periodic. Second, the transformation (2) cannot change the spectrum of a finite- 
dimesional system, so that the deformed chain will have the same spectrum as the 
undeformed one. On an infinite spin chain, we have in principle the possibility to 
use a singular transformation, which will thus change the spectrum. The generators 
of the symmetry algebra J" will transform in the same way as the charges, 

r{x) = s{x)r{o)s-\x) . (3) 
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so the deformation (2) is a morphism of the symmetry algebra^. In particular, if 
the Tindcformcd spin chain has Yangian invariance, this should be also the case 
for the deformed spin chain. The existence of an exact Yangian symmetry im- 
plies that there exists a monodromy matrix satisfying the Yang-Baxter equation 
with a rational i?-matrix. In this case, we can use the algebraic Bethe ansatz 
method to diagonalize the transfer matrix and construct the eigenvectors and their 
scalar products. An alternative way to proceed is to define the long range spin 
chains non-perturbatively in terms of the deformation parameter, as compared to 
the perturbative definition implied by (1). One example of chain defined non- 
perturbatively is the Inozemtsev model [5]. It is known that this model, defined 
with periodic boundary conditions, is diagonalizable by (asymptotic) Bethe ansatz 
in the long chain limit. In the finite size limit the solution is more complicated, 
due to wrapping interactions. The absence of an exact solution in terms of Bethe 
ansatz is a sign that the Yangian symmetry algebra is not preserved by the finite 
size periodic chain. In [1] we have defined a model which coincides in the bulk, in 
perturbation, with the Inozemtsev model but it differs from it by some defect term. 
The defect interaction is long range, with the range growing with the perturbation 
order, and effectively suppresses the wrapping interactions. As a consequence, the 
defect Inozemtsev spin chain preserves exactly the Yangian algebra and is diagonal- 
izable by Bethe ansatz at any value of the deformation parameter and the length. 
The eigenstates and the scalar products are therefore computable by the usual 
algebraic Bethe ansatz procedure, as shown in [1]. In this note, we work out the 
first non-trivial order in the defect term of the Hamiltonian and the corresponding 
Bethe ansatz equations up to next order in the deformation parameter. For long 
spin chains and low perturbation orders, the contribution of the defect is negligible 
and we obtain back the bulk Inozemtsev quantities. We compare with the results 
obtained from the boost deformations of the finite periodic chain. 

The model. In [1] we have studied a long range Hamiltonian which can be 
constructed from the Dunkl operators [6, 7] 

L i-l L i-1 

j=i+l j=l j=l 

where Qij = Zi/(zi — zj) and Zj = e^^'^ with k a real number. The operators 
permute the coordinates KijZj = ZiKy. Under the coordinate permutations, the 
Dunkl operators obey^ the relations of a degenerate affine Hecke algebra [8] 

[di, dj] = , Ki^i+idi — di+iKi^i+i = 1 , 

[Ki,i+i,4]=0 if kj^i,i + l. (5) 

The model is defined via the monodromy matrix 

L 

Ta{u) = 7v{fa{u)) , f„(u) = l[{u- i/2 -id,+ iPja) , (6) 

J=l 

^We thank I. Rostov for this observation. 

^In the defining relations of the degenerate affine Hecke algebra we use no index for the Dunkl 
operator, since they do not depend on the specific representation. 
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where the projection tt transform coordinate permutations Kij into spin permuta- 
tions Pij at the right of an expression, 

^{...Kij) = 7r{...)Pij. (7) 

Using the defining relations (5) of the degenerate affine Hecke algebra, one can 
prove that the Bcrnard-Gaudin-Haldane-Pasquier (BGHP) [9] projection defined 
above is a morphism of the Yangian algebra 

7v{fa{u)fa>{v)) = 7v{fa{u))7V {fa> (v)) , (8) 

and that Ta{u) satisfies the Yang-Baxter equation with the rational R matrix 
R = u-\- iP. The Hamiltonian associated to the model defined above is integrable, 
and when L is large it coincides, in the bulk, with the Inozemtsev Hamiltonian [5] 

When K is large, the interaction falls off rapidly and only the nearest neighbors 
interact. The k ^ oo limit coincides thus with the Heisenberg model. We can use 
K as a tunable parameter to deform the XXX spin chain into a long range spin 
chain. This was done in [10] to reproduce the first three non-trivial orders in the 
perturbative expansion of the dilatation of the J\f = 4 SYM theory in the su{2) 
sector. The deformation parameter is given by the identification 

where IQir'^g^ = A, the 't Hooft coupling constant of the gauge theory. Since 
the sum in (4) is finite and the coefficients Qij are not periodic in i i + L, 
the Hamiltonian obtained from the above Dunkl operators is not translationally 
invariant. Instead, it can be seen as a closed chain with a defect around the sites 
i = 1 and i = L. More general long range spin chain could be built in the same 
manner if we can find other representations of the Dunkl operators. In general, 
long range spin chains can be seen as reductions of short-range interaction systems 
by freezing some degrees of freedom. One particular example is the derivation 
of the BDS spin chain [11] from the Hubbard model at half filling [12]. In that 
example, one obtains a spin system with inhomogeneities 6i = 2g sin qi where 
Qi are the momenta of the underlying fermions. The Dunkl operators can be 
thus seen as some dynamical impurities. On the technical side, the interest in 
working with the operators (4) is that they allow to fully exploit the algebraic 
Bethe ansatz formalism to construct the (cigcn) vectors and their inner products. 
The construction of the eigenvectors and of the conserved quantities for the long- 
range model (6) is facilitated by the link [1] with the inhomogeneous XXX model 
with monodromy matrix 

L 

To{u; 0) = '[[{u- 9j - i/2 + iPja) . (11) 

3=1 

The relation between the two objects is realized by the operator defined by 
T{u) = %To{u;9)\^^^ (12) 



n jl<...<jn kl,...,kn 



=0 



Jl Jn 
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where we have used the notation^ 



(13) 



The morphism property of the BGHP projection (8) gets transferred to the oper- 
ator ^0, 



Ta{u)Ta'{v) = % To,a{u;e) ^^^^^g T^^a'^U-Q) 



%[To,a{u;e)To,a'{v\e)] ^ .(14) 



Upon acting on the pseudo- vacuum state jfi) = | tt • • • t)i the BGHP morphism 
becomes a differential operator, close to the theta quasi-morphism^ defined in 
[15, 16], 



%[...\{u-e) m=ve[...\{u-e) 

U = \J 



(15) 



where the dots stand for any product of elements of the monodromy matrix, and 

' /c X -|- . . . kn 



c: 



ki,...,kn afei 



n ji<...<j„ki, 



Jn 



(16) 



(17) 



The coefficients C'jj^'' computed in principle from the knowledge of the 

Dunkl operators. In [1], we have computed their values in the bulk up to order 
in the perturbative expansion (10), which are translationally invariant. If we take 
into account the effect of the non-periodicity of the definition (4), we have 



L-l 



(18) 



From the above identification we deduced that the (eigen) vectors and their scalar 
products for the long range model can be straightforwardly computed from the 

corresponding quantities in the inhomogcneous model. Denoting as usual with 
A{u), B{u), C{u), D[u) the elements of the matrix T{u), the eigenvectors of the 
transfer matrix, TrT(ti) = A{u) + D{u), can be constructed as 



\{u})g = B{ui) . . . B{uM)\n) = Ve\{u; e})e=o , 



(19) 



where Kuj^}) = Bq{ui;9) . . . Bo{um',0)\^)- Since T{u) obeys the Yang-Baxter 
equation with the rational R matrix, R{u) = n + iP, the algebra of the matrix 
elements is the same as for the usual XXX model, 



A{v)B{u) = ^^—^B{u)A{v) —B{v)A{u) , 



u — V 



u — V 
i 



D{v)B{u) = ^ ' Biu)D{v) + —B{v)D{u) . 



u — V 



u — V 



(20) 
(21) 



^One should keep in mind that in general d, • dj = Tz{di)Tz{dj) ^ d^d^ = ixididj). 
^We call the action of a quasi-morphism because it contains cross terms which spoil the morphism 
property. 
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Therefore, the vectors are eigenvectors of the transfer matrix TrT(u) with 

eigenvalue 



M 



t{u) = a{u) '^^Z...'^ +di^) '^^Z^.^^ > with Q{u) =l[{u - m) (22) 



Q{u) 



Q{u) 



i=l 



provided that the rapidities {u] are satisfying the Bethe ansatz equations 



a{Uj) Uj —Uk + t 



(23) 



Above, a{u) and d{u) the are eigenvalues of A{u) and -D('u) on the pseudo vacuum 



a(u) = /(n + V2)^ , d{u) = f{u - z/2)^ , 



/(k)^|0) = [](n - 0,) |Q)|e=o = 



JJ(M-idj) 



(24) 

\^) . (25) 



To determine the function f{u) we only need the value of the projection of the 
symmetric sums, 



with Cn determining the expansion 



Aln/(n) = i^% 
du u ^ 

n>0 



(26) 



(27) 



By explicit computation using the expression (4) and the perturbative expansion 
(10) we get C2n+i = and Co = 1, C2 = 2^2(1 - i/L) - 2^4 /L + 0(/) and 
C4 = 6(7^(1 — 5/3L) + 0(5^). The part gives the bulk contribution, which was 
already evaluated in [1]. Denoting with x{u) the bulk part, which up to order 
is the same for Inozemtsev and BDS spin chain, we find that 



/« = x« 1 + 



(28) 



It is remarkable that the correction docs not depend on L. Up to the specified 
order, the Bethe ansatz equation can be written as 



X 



+ 



7;.- — 111. — 7' 



/ 3i^k 
with the phase shift (t)o{u) given by 



Uj — Uk — I 



(29) 



1 



1 \ ^g^^i f 1 



{u+f {u-fj 2 V(^+)^ (^") 



+ 0(/), (30) 



where we used the notation = u ^ i/2 and x^ = a;(u^). The phase shift 
4)o{u) in the Bethe ansatz equations (29) can be interpreted as a scattering on a 
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defect situated near the closing point of the chain, k = L. This interpretation 
will become clearer when wc shall derive the explicit form of the Hamiltonian, see 
below. A similar phase shift appeared in the treatment of long-range deformations 
of boundary spin chains, [4] and it can be traced back to the odd conserved charges 
of the chain [2]. Since the strength of the defect does not depend on L (at least 
to the given order), on large chains the effect of the defect will be sublcading, and 
at the leading order in L we retrieve the bulk Bcthc ansatz equations. Let us note 
that equations (29) are exact for any L, provided that we have determined (pQ. For 
the Inozemtsev model, where we know the exact expression of the Dunkl operator, 
this can be done in principle, but we are not doing it here beyond three loop. 

Once the algebra of the matrix elements of the monodromy matrix is known, 
the scalar products of an eigenstate with an arbitrary off-shell state (i.e. with 
rapidities not satisfying Bethe ansatz equations) can be computed by the Slavnov 
method [13], and it reduces to the determinant 



M 



i{M\{u})g = = Yl a{vj) d{uj 



with 



Q,{u, v) 
t{u) = 



t[u 
i 



detofc 



-Vk+i 



u{u + i) ' 



t{v 



(31) 



L M 

n 



Uj + I 



u — It,- 



When {u} = {v} this formula gives the norm of the states in the algebraic Bethe 
ansatz normalization. 

The conserved quantities. On an infinite lattice, the Inozemtsev Hamilto- 
nian (9) can be built from the quantum determinant of the monodromy matrix, 
as explained in [14]. On a finite size lattice the Hamiltonian does not commute 
anymore with the Yangian, so it should not be obtained from the quantum de- 
terminant, but rather from the transfer matrix, by taking higher derivatives^ of 
Inr(n) at = i/2. When = we have the usual XXX conserved Hamiltonians 



•n-l 



^OrilnTo(u)|„=,/2. 



or, explicitly, for the first few charges 

Hi = ln[/o , Uo = TraPai ■ ■ ■ PaL = Pl,l-i ...Pi2 



(32) 



(33) 



H2 = Uq^ ^TraPal...Pak... PaL = ^kk+l , 



(34) 



k j 
where the check means the corresponding factor is absent. The next charges are 



similar strategy was taken in [16]. 
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given by the higher order derivatives 

L 

H3 = y^[-Pfcfc+i;-Pfc+i,fc+2] — L , 
k=l 

L L L 

H4 = ^[-Pfcfc+l, [Pk+l,k+2, Pk+2,k+3]] - ^ Pk,k+2 + 2 ^ Pk,k+1 ■ 
k=l 1=1 i=l 

The eigenvalues of the transfer matrix are given by 

, . , Q{u-i) Q{u + i) 

to{u) = ao{u) h do(u) 



Q{u) 



Q{u) 



(35) 



(36) 



with ao('u) = {u + , do{u) = {u - , (37) 

so that the eigenvalues of the nth Hamiltonian (except for n = 1) are given by 



M 



n-l 



n— 1 



(38) 



where we used the notation = u ± z/2. Now we will turn to computing the 
conserved quantities of the deformed model (6). When the c/j's are present, we 
have 



[/ = TT (TV„ {Pal -di)... {PaL - dL)) , 



(39) 



H2 = idu Inr(ix) = TT C/-^ V TV„ (P„i - di) . . . {Paj - dj) . . . {PaL " C^l) 

u=i 2 \ ' 

\ J 

Computing these operators by evaluating the projections explicitly might be in- 
volved . We are interested here in understanding the structure of the deformation 
in the perturbative parameter from the equation (10). Expanding to the second 
order in the Dunkl operators we obtain an expression which is similar to the Hamil- 
tonian with impurities [3], since before projection the Dunkls behave as c-numbers. 

[/ = f/o j 1 - ^ Pkk+idk + Pkk+iPjj+idkdj + 0{d^) j , (40) 

\ k k<j J 

-f^2 = ^ Pkk+1 + ^ dfc - ^[Pfc_i^fc, Pkk+l] dfc 
k k k 

+ y~l(-Pfcfc+i + Pk,k-i - Pkk+2) d| + ^[[Pfc-ife, Pkk+l], Pk+ik+2] dfedfe+i -I- O(d^) 

k k 

The result of the projection of the Dunkl operators can be computed by hand for 
the terms, and using Mathematica for the y** terms. Up to terms of order g'^ in 
the bulk we have 



dfe — 9^{Pkk-i - Pkk+i), d| — -2g^, dfcdfc+i — 5^ 



(41) 
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and at the boundary we have 

di = -g^Pi2 = dfi'^ - q^Plu dL = g^PL-i,L = d^"''^ + g^PLi, 

d? = di = -g^ = d^''-!'^ + g^, didi = = (didi)^"!'^ - g^ . (42) 

Replacing these values in (40) we get 




H2 = -ff 2 + 9^ ^[-Pfefc-l) [Pkk+l, Pk+l,k+2]] — ^ Pkk+1 
k=l k=l 

+ 2g^PLi-g^[PL-i,L[PLi,Pi2]] + 0{g''). (43) 

If we now combine with higher XXX Hamiltonians, necessary to get rid of the 
multispin bulk interaction, we get 

InC/def = lnU-g^Hl-{2L-l)g^ = lnUo-g^[PLi,Pi2] + 0{g^), (44) 

L L 

H2,dei = H2-g''H2 = {l-4g^)J2Pkk+i+g^Y.^kk+2 

k=l k=l 

+ 2g''PLi-g''[PL-i,L[PLi,Pi2]] + 0{g^) . (45) 

The operator (44) is the momentum operator, while the Hamiltonian (45) coincides 
with the bulk two-loop dilatation operator plus a defect term. The same bulk 
expressions can be obtained by deformation of the spin chain using the procedure 
from equation (1), 

L-l 

X(5') = B[Hl] + 0{g^) = i m,k+i^Pk+i,k+2] + 0(5') , (46) 

k=0 

where B[H^] is called the boost operator of H^. As it is defined above, the boost 
operator also generates a defect term, but which is L times the defect term in (44) 
and (45). This can be understood from the fact that the deformation (1) is in fact a 
similarity transformation and it should not change the spectrum/Bethe equations 
on a finite chain. Let us now characterize the spectrum of the conserved quantities 
of the defect Hamiltonians determined above. For the deformations we consider, 
d{i/2) = a;(0)^, so the conserved quantities, for long spin chains, have the same 
functional form in terms of rapidities as in formula (38) of the XXX spin chain, up 
to some higher order corrections that we can neglect here.^. This does not mean 
that the spectrum of the conserved charges is the same; the quantization conditions 
implied by the Bethe equations are different. Moreover, we have redefined the 
conserved quantities in (44) and (45) by taking linear combinations, so we get 

E...^E2-g^E.^Oig^) ^ C - (^-^ + ^ - - - ^ j + C.(.^) 
^In the BDS spin chain the corrections are of the order . 
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with C = iL{f{i) + 2g^ f"[i)). In this way we obtain a perturbative expression 
which matches the charges for the BDS spin chain [11]. Prom (44) and the prop- 
erties of the shift operator Uq we conclude that 

(C/def)'' = 1 - 9\H^ + L) + 0{g^) . (48) 
Let us now look at the one magnon states with wave function given by 

L 

^p = Y.{e'^' + g'ck)\k). (49) 
fe=i 

where we have denoted with \k) = a~^\^) . The coefficients are zero in the bulk. 
Their value near the defect can be computed in the simplest way by imposing that 
*p is an eigenvector of U^^i with eigenvalue e^^. We obtain 

CL = e^'P° - e'P° , Cl = e-'P° - 1 . (50) 

From the previous formula, we obtain, at two loops, 

pL = 2nn + ig' - ^) + 0(/) (51) 

so that the quantization condition for the momentum p is 

pL - 16g^ cosp/2 sin^p/2 = 2™ , (52) 

1%^ , o , , 27rra 

or P = Po -\ p — cospo/2sm ]9o/2 , with po = —— . 

At the next loop order, taking into account the expression (30) of the defect phase 
shift one has 



, I65 . 3P0 Po 

P = Po -\ sin — cos — 

^ ^ L 2 2 



1— 5^ (4— 10cospo+5cos 2po + ^(cos 2po—cospo 



The vector diagonalizes also the defect Hamiltonian (45) with eigenvalue given 
by (47). Let us emphasize that, although the functional form of the eigenvalues 
(47) are the same as in the translational invariant case, the actual eigenvalues are 
modified due to the change of the Bethe equation (quantization condition in the 
case of the one- magnon states), 

E2MP) = -4sin2p/2 + 16g2sin4p/2 + 0(54) (53) 

= -4sin2j9o/2 + 16c/^sin^j9o/2 ^ cos^j9o/2sin^po/2 + 0(9"^) • 

We have checked that the two-loop eigenfunctions can be obtained alternatively by 
applying the operator Vq to the eigenvectors of the inhomogeneous XXX model, 
cf. equation (19). 

Conclusion We have studied the conserved quantities of a long-range defor- 
mation of the XXX model that is defined globally, i. e not only perturbatively. This 

model coincides in the bulk with the Inozemtsev model [5], which in turns coincides, 
up to three- loop order in the deformation parameter, with the BDS model [11] and 
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with the dilatation operator on the A/" = 4 SYM theory. The difference with the 
the Inozcmtscv model is given by a defect term situated at the closing point of 
the chain. The exact spectrum and eigenvectors of the resulting Hamiltonian can 
be built by algebraic Bethe ansatz, for any length of the spin chain. On a chain 
of length L they differ from the corresponding quantities of the periodic chain by 
terms of order 1/L. This justifies their use in [1] to extend the semiclassical results 
of [17] to two-loop three-point correlation functions in the su{2) sector of TV = 4 
SYM [18, 19, 20]. 
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